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PAPER  NO.  36 


A  MATHEMATICAL  MODEL  FOR  THE  BEHAVIOR  OF  THE  BRAIN 
WHEN  THE  HUMAN  HEAD  IS  SUBJECTED  TO  IMPULSIVE  LOADS 


V 


by 

All  E.  Engin,  Ph.D. 
and 

Verne  L.  Roberts,  Ph.D. 

The  University  of  Michigan 

ABSTRACT 

The  subject  matter  of  this  paper  Is  concerned  with  the  theoretical  deter¬ 
mination  of  the  behavior  of  the  brain  when  the  human  head  Is  subjected  to 
external  Impulsive  loads.  The  mathematical  analysis  Is  made  for  the  axlsymmetrlc 
response  of  an  Inviscld  compressible  fluid  loaded  Impulsively  by  Its  elastic 
spherical  shell  container.  The  motion  of  the  fluid  Is  assumed  to  be  governed 
by  the  linear  wave  equation.  The  spherical  shell  equations  Include  both 
membrane  and  bending  effects  In  axlsymmetrlc  torsionless  motion.  In  the 
analysis  first  the  solution  for  an  Initial  value  problem  is  obralned;  later 
the  transient  response  of  the  fluid  for  an  arbitrary  velocity  Input  of  the  shell 
Is  constructed  by  means  of  convolution  Integral.  For  the  numerical  results,  a 
characteristic  time  Is  defined  and  the  excess  pressure  distribution  In  the 
fluid  Is  evaluated  for  various  deceleration  times  comparable  with  this 
characteristic  time.  A  description  of  some  of  the  salient  features  of  the  excess 
pressure  distribution  Is  also  given  In  view  of  the  elastic  and  rigid  shell 
boundary  of  the  fluid.  Since  the  problem  is  considered  from  a  theoretical 
point  of  view  to  determine  some  of  the  causes  of  the  brain  damage  when  the  human 
head  Is  subjected  to  Impulsive  loads,  in  the  numerical  computations  the  data 
Is  chcsen  to  be  suitable  to  the  physical  properties  of  brain  and  skull. 
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LIST  OF  SYMBOLS 

E  Young's  modulus 

Pn(cos  $)  Legendre  polynomials  of  the  first  kind 

P^(cos  $)  Associated  Legendre  polynomials  of  the  first  kind  and  first  order 
V  Speed  of  the  fluid-filled  shell  along  an  axis  passing  through  the 
poles 

♦  Velocity  potential  ft**  the  fluid 

♦j  Nondimensional  velocity  potential  for  the  fluid,  i>/acs 
n  Nondimensional  frequency,  «a/c 
a  Radius  of  spherical  shell 

ao’an  Coefficients  of  Legendre  polynomial  expansion  of  t 
bn  Coefficients  of  Legendre  polynomial  expansion  of  y 
co,cn  efficients  of  velocity  potential 
c  Compresslonal  wave  speed  in  the  fluid 
c$  Apparent  wave  speed  in  the  shell,  [E/ps(1-v2)]j/2 
f  Shell-fluid  parameter,  Poa/p$h 
h  Shell  thickness 

j„(z)  Spherical  Bessel  function,  {if/2z)'l/20(z) 
n  n+1/2 

k  Wave  number,  w/c 

p  Excess  pressure 

pt  Nondimensional  excess  pressure,  p/pQc2 
r,e,$  Spherical  coordinates 
Tj  Nondimensional  radius,  r/a 
s  Speed  ratio,  c/c$ 
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t  Time 


*1 

tc  Characteristic  time,  ^r 
tQ  Deceleration  time 

u  Meridional  displacement  of  the  shell  mid-surface  with  respect  to 

J 

geometric  center  of  the  shell 

I 

w  Radial  displacement  of  the  sh*ll  mid-surface  with  respect  to  geometric 
center  of  the  shell 
o2  Thickness  parameter,  h2/12a2 
C  Nondimensioiial  radial  displacement,  w/a 
i>  Nondlmenslonal  meridional  displacement,  u/a 
An  n(n+l) 
v  Poisson's  ratio 

p  ,ps  Mass  density  of  fluid  and  shell  respectively 

t  Nondimensional  time,  cst/a 

t  Nondlmenslonal  characteristic  time,  2/s 
c 

t0  Nondimensional  deceleration  time,  c$t0/a 
u  Angular  frequency 

I 
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INTRODUCTION 


This  Investigation  alms  at  the  following  two  considerations.  First,  the 
subject  matter  Is  a  point  of  Interest  In  theoretical  mechanics  due  to  the  fluid- 
solid  Interaction  nature  of  the  problem.  Second,  from  the  application  point  of 
view,  the  Impulsive  response  of  the  fluid  when  the  enclosing  elastic  shell  Is 
suddenly  subjected  to  a  change  In  Its  velocity  can  be  taken  as  a  simple  but  Im¬ 
proved  theoretical  model  to  determine  the  formation  of  brain  damage  when  no  local 
contacts  are  made  on  the  human  head. 

The  previous  studies  In  the  area  have  been  either  Investigations  Involving 
shells  In  contact  externally  and/or  Internally  with  fluids,  especially  by  researchers 
In  the  field  of  acoustics,  or  analyses  of  various  head  Injury  models.  While  studies 
on  these  two  categories  are  numerous,  only  a  few  representative  ones  will  be  men¬ 
tioned  here.  Junger1  Investigated  the  effect  of  fluid  on  the  natural  frequencies 
of  cylindrical  and  spherical  shells  freely  suspended  In  a  compressible  fluid  medium. 

Free  and  forced  oscillations  of  Infinitely  long  cylindrical  shells  surrounded  by 
water  were  studied  by  Greenspon2  who  treated  unpressurized  shells  by  exact  elasticity 
theory  and  cylindrical  shells  with  fluid  by  approximate  shell  theory.  Utilizing 
linear  shell  theory,  which  Includes  both  membrane  and  bending  effects,  Engin  and  Liu3re- 
cently  obtained  the  frequency  equation  and  corresponding  frequency  spectrum  of 
fluid-filled  spherical  shells  for  the  axisymmetric  and  nontorsional  motion. 

In  the  literature,  the  rigorous  mathematical  treatments  of  the  physical  theory 
of  the  formation  of  brain  damage  was  first  Introduced  by  Anzelius4  and  Giittinger5 
with  their  analyses  of  the  Impulsive  response  of  an  inviscid  fluid  contained  In  a 
rigid  closed  spherical  shell  (or  container).  Their  formulations  are  essentially 
Identical  and  Involve  an  axisymmetric  solution  of  the  acoustic  wave  equation  in 
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spherical  coordinates.  In  the  papers  of  both  authors  the  eigen  values  of  the 
problem  are  determined  by  requiring  the  radial  component  of  the  fluid  velocity 
to  vanish  at  the  Interior  surface  of  the  rigid  spherical  shell  surrounding  the 
fluid.  Hayashl6  treated  a  one-dimensional  version  of  the  Anzellus-Giittinger 
model.  His  model  consists  of  a  rigid  and  massless  vessel  containing  Inviscid 
fluid.  The  vessel, which  Is  attached  to  a  linear  spring,  is  subjected  to  impacts 
with  a  stationary  wall.  Approximate  solutions  were  obtained  for  the  limiting 
cases  of  soft  and  very  hard  Impacts.  Although  this  simple  model  has  the  advan¬ 
tage  of  being  easy  to  interpret,  It  has  the  similar  shortcomings  of  the  Anzelius- 
Giittlnger  model,  l.e.  (a)  due  to  rigidity  and  geometric  assumption  there  is  no 
way  to  determine  the  possible  locations  of  skull  fracture  and  (b)  the  effects  of 
skull  deformation  on  the  Intracranial  pressure  distribution  can  not  be  determined. 
Recently,  Engin7  removed  some  major  restrictions  of  previous  head  injury  models 
by  obtaining  analytical  and  numerical  solutions  for  the  dynamic  response  of  a  fluid- 
filled  elastic  spherical  shell.  The  loading  pattern  for  his  modr l  is  taken  to  be 
local ,  radial ,  ax 1 symmetric  and  impulsive.  Since  the  load  is  applied  as  a  force 
locally  on  one  of  the  poles  of  the  shell,  the  combined  shell  theory  which  includes 
both  membrane  and  bending  effects  of  the  shell  has  been  used  for  the  proper  de¬ 
scription  of  the  wave  propagation  on  the  shell.  The  analysis  utilizes  Laplace 
transform  technique  in  obtaining  the  transient  response  of  the  system.  The  con¬ 
clusions  of  Engin's  paper  include  the  possible  locations  of  brain  damage  and  skull 
Injury  on  the  basis  of  the  numerical  computations.  As  a  problem  in  mechanics, 
the  present  investigation  is  a  generalization  of  the  results  of  Anzelius  and 
Guttinger  by  removing  the  restriction  of  rigidity  of  the  shell  surrounding  the 
fluid.  Our  model  consists  of  an  elastic  spherical  shell  filled  with  inviscid  com¬ 
pressible  fluid.  The  shell  material  and  fluid  arc  considered  to  be  homogeneous  and 
isotropic.  In  the  analysis,  the  fluid-filled  shell  will  be  considered  to  have  a 
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constant  translational  velocity  for  t<o  with  respect  to  an  Inertial  coordinate 
system.  At  tco  the  shell  Is  brought  to  a  sudden  stop,  l.e.  the  fluid  occupying 
the  Interior  space  of  the  shell  Is  subjected  to  a  global  axisymmetric  Impulse 
on  Its  boundary.  The  determination  of  the  pressure  distribution  In  the  fluid 
for  this  kind  of  Impulse  will  help  to  explain  quantitatively  the  location  and 
the  magnitude  of  brain  damage  under  the  conditions  In  which  the  application  of 
local  forces  on  the  skull  is  avcIJed. 
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I.  EQUATIONS  OF  MOTION  AND  THEIR  SOLUTIONS 


4 


The  governing  differential  equations  of  motion  for  a  fluid-filled  spherical 
shell  were  previously  obtained  In  reference  3  by  means  of  Hamilton's  Principle. 
These  three  partial  differential  equations,  which  are  coupled  In  terms  of  the 
meridional  and  radial  displacements,  u,  w  of  the  shell  mid-surface  and  velocity 
potential,  t,  of  the  fluid,  are  given  below  In  nondlmenslonal  form: 

°2[f*£  +  cot*  S*  ‘  "  ||§-  -  cot*  +  (v+cot?-4>)  |J-j+  +  cot*  |J 

-  (v+cot2*)*  +  (1+v)  -0,  (1) 

°2[fff  +  2cot+  "  0+v+cot2*)  fj-  +  (cot2<>-v+2)^cot*  -  —£  -  2cot$  ||f 
♦  (1+v+cot2*)  -  (2-v+cot24>)cot*  |^-j-  (l+v)(||-+  *cot<f>+2c)  -  |^§- 

-  f  *  0,  (2) 


Since  the  deformation  of  a  given  shell  Is  usually  analyzed  terns  of  the  de¬ 
formation  of  Its  mid-surface,  Eqs.  1  and  2  describe  the  motion  of  an  imaginary 
boundary  of  the  fluid  which  is  half  the  shell  thickness  av/ay  from  the  physical 
boundary.  These  shell  equations  Include  both  membrane  and  bending  effects  and 
they  are  given  for  an  axi symmetric  torsionless  motion.  Equation  3  is  the  linear 
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wave  equation  describing  the  motion  of  small  oscillations  of  Invlscld  and  1r- 
rotatlonal  fluid. 

Consider  the  following  series  expansions  for  the  nondimenslonal  radial  and 
tangential  displacements  of  the  shell  midsurface: 

c(*,t)  *  l  an{t)  Pn(cos$),  (4a) 

nco 

and 


■  l  bn(T)  P^cos*),  (4b) 

n=l 

where  Pn(cos<*.)  are  Legendre  polynomials  of  the  first  kind  and  P^(cos«j>)  are  asso¬ 
ciated  Legendre  polynomials  of  the  first  order,  first  kind.  Since  the  second 
solutions  of  the  Legendre  equations  are  singular  at  the  poles  they  are  not  in¬ 
cluded  In  the  expansions  of  c  and  <|>.  The  requirement  of  boundedness  of  solutions 
and  the  linearity  of  Eq.  3  lead  to  its  formal  solution: 

CO 

♦i (r,  ,4>,t)  *  l  cn(nr)  jri(kar,)  Pn(cos$),  (5) 

n=o 

where  Jn(kar,)  is  spherical  Bessel  function,  kew/c  is  the  wave  number  and  to  Is 
the  circular  frequency. 

The  boundary  condition  between  the  fluid  and  shell  can  be  stated  as  the 
continuity  of  normal  velocities  for  all  $  and  i;  that  is, 

isl&ill.  B  (1  (gj 

3t  ar  i 


Substitution  of  Eqs.  4a  and  5  into  Eq.  6  yields  the  following  relationship 


between  an(x)  and  c0(t)  for  each  n: 

/  \  _  1  dan(T) 
cn(t)  1 QjJW  dT - 


(7.) 
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It  can  be  shown  that  substitution  of  Eqs.  4a,  4b,  and  Eq.  5,  together  with  Eq.  7 
Into  Eqs.  1  and  2,  and  the  repeated  utilization  of  the  differential  equations 
satisfied  by  Pn  and  yields  the  following  system  of  equations  for  the  deter¬ 

mination  of  an(T )  and  bn ( t ) : 
for  n«o: 

r  Jft(n)  1  d2a  (t) 

['  ♦  f  aJ^Tj  V2>  ■  °>  <8> 

for  n  i  1: 
d2b  (t ) 

+  [lW(1-v-xn)]an(,)  -  (l-v-xn)(Ha»)bn(T)  a  0,  (9) 

T  Jn(nf|  <l2«n(i)  (  1  f 

['  *  fW^r\  3^ - ltu")xn  +  «*  tx*-xn(l-v)3j  b„(,)  ♦  J2(Hv) 

•  +  o2t»J-xn(l-v)]J  a„(t)  •  0  (10) 

where  and  xncn(n+l).  In  reference  3  a  description  of  some  of  the  salient 

features  of  the  frequency  spectrum  of  a  fluid-filled  spherical  shell  Is  given  In 
detail.  There  It  was  also  shown  that  for  each  mode  number,  n,  the  infinite  num¬ 
ber  of  frequencies  (or  characteristic  roots)  exists.  Thus,  the  solutions  of  Eqs. 

8,  9  and  IQ  can  be  written  In  the  following  form: 

CO 

a-(-r)  =  l  A-_  sin(n  ST+a_m) ,  (11) 

0  om  om  om'  *  ’ 

m=o 


'  l  Anm  si,l(n 


m=o 


nm 


$T+a  ) 

nm' 


and 


(12) 


bn(T) 


l 

m=o 


U  A™  s1  n(n, 


nm 


S5+anm> 


(13) 
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where  A^,.  A^,.  ^  and  am  are  arbitrary  constants;  and  nnm  are  the  roots 
of  the  following  frequency  equations  respectively 


f  40(n) 

1  *  f  53JlnTj 


s2n2  -  2(1+v)  ■  0  for  n»o, 


04) 


and 


1  •  ’ajS]  ■**"  • 

{ 


-  (1+v)^  2(1-v-Xn)(1+a2)  +  Xn[l+v-a20-v-Xn)]  j  “ 


a2(2-Xn)[X2-Xn(l-v)]  -  0 


for  n  a  1.  (15) 

One  of  the  steps  In  derivation  of  Eq.  15  also  yields  the  following  expression 

for  w 


"  s^n’A'fl^u-'xJJ  i  ( 1  +a2 ) 


n*1,2,3,... 

m=0,l ,2». .. 


In  view  of  Eqs.  7,  11  and  12  the  nondlmenslonal  velocity  potential,  4>,  now  can 
be  written  as 

-  j  |  «  jjrr  cos(n™  ST+«nn>>  Pn(ctB*>  <16> 

n*o  m*o  n  nm' 

Next,  let  us  consider  the  following  physical  situation  in  which  a  fluid- 
filled  shell  Is  travelling  with  a  constant  speed,  V,  along  an  axis  passing  through 
the  poles.  Assume  at  t*o  It  Is  brought  to  a  sudden  stop.  Incidentally,  this 
situation  can  be  visualized  as  the  motion  of  the  fluid-filled  shell  in  a  force 
field  which  can  only  affect  the  shell  material;  the  parameters  of  the  force  field 
can  be  adjusted  In  such  a  manner  that  the  desired  motion  of  the  shell  Is  obtained. 
Due  to  Inertia  of  the  fluid  particles,  at  time  t=o  they  will  experience  a  velocity 
relative  to  the  coordinate  system  whose  origin  is  located  at  the  geometric  center 
of  the  shell.  In  other  words,  when  the  shell  surrounding  the  fluid  is  brought  to 
a  sudden  stop  at  t=o,  the  fluid  particles  occupying  the  interior  space  of  the 
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shell  are  unaware  that  the  motion  of  shell  Is  arrested.  Assuming  that  the  Initial 
excess  pressure  distribution  In  the  fluid  Is  zero  we  can  now  write  the  following 
Initial  conditions  on  the  nondlmensional  velocity  potential 


♦.(r,,'f,T»o)5*0«^-  r,  cos* 


and 

111  e  *•  ■  n 
St  *o  w 
t*o 

From  the  second  Initial  condition,  l.e.  tj^o,  we  get  0^*0.  The  first  initial 
condition  Is  used  to  obtain  the  coefficients  Ann>.  In  order  to  do  this  let  us 


(17) 


first  write  Eq.  16  In  the  terminology  of  eigenfunctions  by  defining  e1 

VQt»r'>pn<e°s*>; 

«  •  COS(fl  St) 

♦.  (r.  ,*.t)  -  I  I  sAmem  JJtiO - 

n*o  m*o  n'  nm' 


nm 


(18) 


Let  enm“^n^ftnlr,^Pn^C0S*^  be  the  conJu9ate_e^9enfunct1on  for  cnmi  hav*n9  aPPHed 
the  first  of  the  Initial  conditions  on  «•  In  Eq.  18  and  multiplying  both  sides  of 


the  resulting  equation  by  enm  and  Integrating  It  over  the  fluid  volume  we  get 


/u  ♦.  enm  dv 
v  o  nm 


(19) 


Mn  3IW  Gnn  enm 
n°o  m=o  n  ' 


To  evaluate  the  Integral  on  the  numerator  of  A^  v/e  write  cos*  as  P,(cos*)  and 
see  that  In  view  of  orthogonality  of  Legendre  polynomials  in  the  Interval  [-1,1] 


f¥  *o  enm  dv  "  (CC*  c^ar'^3jn^nlr,^P^C0S*^Pn^C0S*^  s1M*dodr,  =  0 

s 

for  n  t  1.  (20) 

Expressing  the  spherical  Oessel  functions  In  terms  of  fractional  order  Bessel 
functions  and  performing  the  Integration  in  Eq.  20  yields 
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where  =  n|m.  For  the  evaluation  of  the  Integral  on  the  denominator  of  Anm  wi 
use  the  following  properties  of  Bessel  functions  and  Legendre  polynomials: 

(*a2[J£+i<knma>  -  Vi(knmaK+f(knma)]  lcm 


■Wknlr)  Jn+i<knmr>rdr  m] 


.  Vm 


and  from  the  definition  of  the  norm  of  Pn(cos$) 

/*  P*(cos$)  sin$  d$  ■  3— 

Thus,  the  value  of  the  integral  under  consideration  becomes 


•  * 

, i (n  )  -  o  i (n  )  0,3 (n) 

r  T  y  s  -  ,  2na3s  n+2  nrn  n-^  r.m'  “n+2*  nm'J 

'  l  JL  JL  3IW  enm  enm  dv  "  2n+l  J  /2  T n  i  \  i 
v  nso  m-o  ‘V  «nld  [  Jn+I(fW  "  °n+|(nnm) J 

From  Eq.  20  v:e  see  that  Anm*0  except  n=l.  Let  A^sA^,  using  Eq;  21  and. 22  in 


Eq.  19  v/e  obtain  the  final  expression  for  Am: 


where  are  the  roots  of  the  frequency  equation  for  n-l.  Substitution  of  Eq. 
23  into  Eq.  18  after  some  simplification  gives  the  expression  for  the  nondimen- 
sional  velocity  potential  for  the  fluid 


cos(V'> 


The  nondirr.cnsional  excess  pressure,  p(,  is  equal  to 
directly  from  Eq.  24. 


Hi 

3t 


and  it  is  obtainable 


It  Is  Interesting  to  note  that  the  appropriate  limiting  case  of  Eq.  24  agrees 
with  the  result  obtained  by  Guttinger5.  When  we  write  Eq.  24  in  dimensional  form, 
c$,  the  wave  speed  on  the  shell,  disappears  from  the  expression  of  velocity  poten¬ 
tial.  In  the  case  In  which  the  shell  surrounding  the  fluid  becomes  rigid,  s-k>, 
and  Eqs.  14  and  15  degenerate  to  j^{fl)co,  v'hich  Is  easily  shown  to  be  the  frequency 
equation  of  an  ideal  fluid  contained  In  a  rigid  spherical  shell.  When  the  shell 
becomes  rigid  n^s  in  Eq.  24  are  taken  to  be  the  roots  of  j|(n)*o  which  can  also 
be  written  in  terms  of  fractional  order  Bessel  functions  as  2nJ|' (ri)-J-j(n)®o. 

This  last  equation  was  used  by  GUttlnger  to  obtain  the  necessary  natural  fre¬ 
quencies  of  the  fluid  for  the  determination  of  the  impulsive  response. 

If  desired,  the  transient  response  of  the  fluid  for  an  arbitrary  velocity 
input,  V(t),  of  the  surrounding  shell  can  be  examined  by  means  of  Convolution 
Integral.  By  the  principle  of  superposition,  It  can  be  shown  that  velocity  po¬ 
tential  o,,  to  anarbltrary  shell  excitation  V(t)  can  be  expressed  as 

♦,(r,,<J,T)  =  *,(rt,$,T)  +  f1  Mr,,*,T-s)de,  for  t>0.  (25) 

where  *, (r ,  ,<> ,t )  is  given  in  Eq.  24  and  V  is  the  constant  speed  before  the  sudden 
stop  of  the  spherical  shell. 


889 


II.  NUMERICAL  RESULTS  AND  DISCUSSION 


For  the  determination  of  the  numerical  values  of  the  excess  pressure,  p, , 
generated  In  the  fluid  we  use  a  different  form  of  Eq.  24  In  which  the  Bessel 
functions  of  fractional  order  arc  expressed  In  terms  of  trigonometric  functions. 
The  expression,  thus  obtained  for  p,  Is 


2Vc  t3(s1rfB-cosn|n)-Vlnf»m][s1r(nn|r,)-fimr,cos(V,)]s1n(omS,) 

V  "Vo  r.,iv ♦  (,-3i>«)s,n“nin  ♦ 


Here  the  dimensional  excess  pressure,  p,  can  be  obtained  from  pcp,p  cc2.  At 
r,co  the  equation  for  p,  Is  Indeterminate!  however,  application  of  L'HospItal 
rule  to  Eq.  26  once  yields  p,(o,$,t)»o.  That  Is,  In  the  equltorlal  plane  which 
is  perpendicular  to  the  direction  of  Impulse,  pressure  Is  zero  at  all  times. 
Since  the  problem- was  considered  to  serve  as  a  theoretical  model  to  determine 
the  formation  of  brain  damage,  the  numerical  values  arc  obtained  from  the  fol¬ 
lowing  data  suitable  to  the  physical  properties  of  brain  and  skull: 


p$  »  0.0772  lbm/in3 
E  *  2  x  106  lbf/in2 
v  -  0.25 

a  B  3  In  (27) 

h  *  0.15  in 
P#  «  0.0362  lbm/in3 
c  B  57100  In/scc. 

The  assumption* leading  to  the  above  data  are  discussed  in  detail  In  references 
7  and  B. 
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Since  the  sudden  stop  of  a  fluid-filled  shell  demands  an  Infinite  deceler' 
ation,  It  Is  more  reasonable  to  consider  the  case  where  stopping  occurs  during 
a  finite  time  tQ  (Nondlmenslonal  value,  t^-j®-).  More  specifically  let  us  as¬ 
sume  that  the  shell  has  the  following  velocity  form 


V(t)  - 


r  V  (tSo) 

<  Vt/t0  (o<t<t0) 


l  o  (tU0)  . 

I’.e.  It  Is  brought  to  a  stop  with  a  constant  deceleration  during  time  t0>  For 
the  numerical  results  the  magnitude  of  V  Is  taken  to  be  528  In/sec.  As 
regards  the  response  of  the  system,  the  above  situation  Is  Identical  to  that 
in  which  the  shell  and  Its  content  are  set  Into  motion  from  rest  with  the  same 
magnitude  of  acceleration.  Thus,  In  view  of  Eq.  25  the  expression  for  the 
excess  pressure  can  be  obtained  directly  from.  Eq.  26 


(28) 


P,(«\.*,t)  *  l  g^.r,)  {[u(T)-u(x-T0)3s1n2 
m=o 

(29) 

♦  u(t-tq)  sin  -^-9-  s1nnws(t~|)| 


where 

„  u  ,  ,  .  4V_  [3(Hn^-cos^)-Vl%][s1n(v,)-V|COi(V,)] 

V*’  1  ct  r  /(L[iii+n-3(U)s1nzn +n„s1na „cos(iJT  cos*> 

so  i  mm  hi  mm  m  m 

and  u(t) ,  u(t-TQ)  are  the  unit  step  functions.  We  note  that  the  limiting  case  of 
Eq.  29  in  fact  gives  Eq.  26  when  t0>o  .  The  series  in  Eq.  29  exhibits  much 
better  convergence  than  the  one  given  In  Eq.  26. 
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let  us  define  a  characteristic  time,  t,.(or  'C„  ■  ,  3.615  for  the  data 

given  In  Eq.  27)  to  be  the  time  required  for  a  wave  in  the  fluid  to  travel  from 
one  pole  to  the  other.  In  all  the  figures  the  deceleration  duration  Is  considered 
with  respect  to  this  characteristic  time.  Figure  1  shows  the  plot  of  the  excess 
pressure  at  the  pole  U»o)  for  two  different  values  of  deceleration  duration.  The 
reduction  of  the  pressure  with  doubling  cf  the  deceleration  time  Is  readily  apparent. 
Figures  2  through  8  present  the  ®xcess  pressure  distribution  along  the  polar  axis 
for  various  values  of  nondimenslonal  radius  r,  and  nondlmonslonal  time  i.  From 
Eq.  29  we  see  that  p,(-r,,$,T)  B  -p^r^M);  hence,  the  points  located  symmetri¬ 
cally  with  respect  to  the  equltorial  plane  which  Is  perpendicular  to  the  polar 
axis  always  experience  the  excess  pressure  of  opposite  sign.  That  Is,  for  such 
a  pair  of  points  If  one  Is  In  a  state  of  compression,  the  other  one  will  be  In 
a  state  of  dilatation.  From  Figs.  2  through  8  one  can  find  the  magnitude  of  the 
excess  pressure  along  any  ray  extending  from  the  origin  by  simply  multiplying  the 
value  with  cos$.  In  other  words,  the  waves  generated  on  the  Inner  surface  of  the 
shell  have  amplitude  factors  which  attenuate  according  to  the  polar  angle  When 
the  deceleration  time  t0«tc,  sharp  extremes  In  the  excess  pressure  form  for  t>tq. 
This  fact  Is  Illustrated  In  Fig.  2.  A  further  Interesting  observation  Is  the 
comparison  of  the  pressure  distributions  In  the  fluid  contained  in  the  rigid  and 
the  elastic  shells.  In  Figs.  3  through  8  pressure  distributions  are  plotted  for 
both  clastic  and  rigid  cases.  From  these  plots  we  see  that  the  amplitude  of  the 
excess  pressure  v/ave  In  the  fluid  contained  in  an  elastic  shell  is  considerably 
loss  than  that  of  the  excess  pressure  wove  in  the  fluid  contained  In  a  rigid  shell. 
Local  variations  of  the  wave  form  In  the  clastic  case  are  also  noticeable.  When 
we  compare  Figs.  5  and  7  and  Figs.  6  and  8  we  notice  that  the  wave  forms  are 
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essentially  the  same  and  the  amplitudes  of  pressure  for  the  large  value  of  t0 
•re  much  less  than  those  corresponding  to  a  small  value  of  tq  (compare  the  scales 
In  the  figures).  In  Figs.  7  and  8  t  «100tr  corresponds  to  502  In/sec2  (or  130  g's 
for  a  deceleration  time  of  0.0105  sec.)  In  this  case  the  maximum  negative  excess 
pressure -generated  In  the  fluid  Is  about  48  psl. 

In  conclusion,  If  we  seek  the  brain  damage  to  occur  at  the  points  of  rare¬ 
faction  of  the  fluid,  we  find  this  situation  arises  In  maximum  magnitudes  at  the 
poles  and  these  locations  are  quite  significant  for  the  analytical  confirmation 
of  the  cavitation  theory  of  brain  damage. 
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Nondimensional  excess  pressure,  pt,  vs.  nondimensional  time,  x,  at  $=o 
(Elastic  shell  boundary). 


10  x60 


Nondimensional  excess  pressure,  p,,  vs.  nondimensional  time,  x,  and 
nondimensional  radius,  r.;  $=o  and  x  =.01x  ,  (Elastic  shell  boundary) 


OSxp..(H 
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Nondimcnsicnal  excess  pressure,  p|t  vs.  nondimensional  time,  t,  and 
nondimensional  radius,  r  ;  *=o  and  t  =0.5x  ,  (Rigid  shell  boundary) 


10"Jx20 


Nondimensional  excess  pressure,  p,,  vs.  nondimensional  time,  t,  and 
nondimensional  radius,  r  ;  4*0  and  t  *0.5x  ,  (Elastic  shell  boundary) 


10‘6x150 


100 
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Pig.  7  Nondlmenslonal  excess  pressure,  p  ,  vs.  nondlmenslonal  time*  x,  and 

1  t 

nondlmenslonal  radius,  r, ;  *«o  and  to»10Otc,  (Rigid  shell  boundary) 
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Nondimensional  excess  pressure,  pt ,  vs.  nondiciensional  time,  t,  and 
nondimensional  radius,  r,i  **o  and  t  “lOOr  ,  (Elastic  shell  boundary) 
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